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Well-posed boundary integral equation formulations and Nystrom 
discretizations for the solution of Helmholtz transmission problems 
in two-dimensional Lipschitz domains 

Victor Dominguez* * Mark Lyonj Catalin Turc^ 


Abstract 

We present a comparison between the performance of solvers based on Nystrom discretiza¬ 
tions of several well-posed boundary integral equation formulations of Helmholtz transmission 
problems in two-dimensional Lipschitz domains. Specifically, we focus on the following four 
classes of boundary integral formulations of Helmholtz transmission problems (1) the classical 
first kind integral equations for transmission problems |13j . (2) the classical second kind integral 
equations for transmission problems [25], (3) the single integral equation formulations [21j . and 
(4) certain direct counterparts of recently introduced Generalized Combined Source Integral 
Equations HIS]. The former two formulations were the only formulations whose well-posedness 
in Lipschitz domains was rigorously established m We establish the well-posedness of 
the latter two formulations in appropriate functional spaces of boundary traces of solutions of 
transmission Helmholtz problems in Lipschitz domains. We give ample numerical evidence that 
Nystrom solvers based on formulations (3) and (4) are computationally more advantageous than 
solvers based on the classical formulations (1) and (2), especially in the case of high-contrast 
transmission problems at high frequencies. 

Keywords: transmission problems, integral equations, Lipschitz domains, regularizing op¬ 
erators, Nystrom method, graded meshes. 
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1 Introduction 

A wide variety of well-posed boundary integral equations for the solution of Helmholtz transmission 
problems has been proposed in the literature, at least in the case when the interfaces of material 
discontinuity are regular enough. Most of these formulations are derived from representations of 
the fields in each region filled by a homogeneous material by suitable combinations of single and 
double layer potentials. The enforcement of the continuity of solutions and their normal derivatives 
across interfaces of material discontinuity leads to Combined Field Integral Equations (CFIE) 
of transmission scattering problems. Some of these integral formulations involve two unknowns 
na EoiEaEaEi], three or more unknowns 123 [33], while others involve one unknown per each 
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interface of material discontinuity [2T]. It is also possible to formulate transmission problems in 
terms of both interior and exterior traces-Multi-Trace Formulations (MTF), that is using four 
unknowns per each interface of material discontinuity mm- More general boundary problems, 
which include not only transmission conditions but mixed Dirichet-Neumann conditions in complex 
geometrical configurations as well, have been studied in [38]. 

In the technologically important case of transmission problems at high-frequencies, the numeri¬ 
cal solutions of boundary integral equation formulations typically rely on Krylov subspace iterative 
methods. As in the case of impenetrable scattering problems |6], the classical boundary integral 
equations of transmission problems are not particularly well suited for Krylov iterative solutions 
of transmission problems at high-frequencies. We have demonstrated recently that a novel class 
of boundary integral equations referred to as Generalized Combined Source Integral Equations 
(GCSIE) [H [5| is a more favorable alternative for smooth transmission problems that involve high- 
contrast configurations at high-frequencies. The main scope of this paper is to investigate to what 
extent the aforementioned claim is valid in the case of high-frequency, high-contrast Helmholtz 
transmission problems when the interface of material continuity is a Lipschitz curve. 

An important question related to boundary integral equation (BIE) formulations of linear, 
constant-coefficient PDEs is whether the BIE are well-posed. This issue is typically settled via 
Eredhohn arguments whose flavor differ significantly from the case of regular boundaries to the 
case of Lipschitz boundaries. The case of smooth boundaries is extremely well understood and 
researched, as one can take full advantage of the increased smoothing properties of double layer 
operators that guarantee compactness properties needed in the Fredholm theory |2j . In addition, a 
very general methodology based on coercive approximations of Dirichlet-to-Neumann operators can 
deliver optimally conditioned boundary integral formulations (GCSIE) for transmission Helmholtz 
problems |5]; the aforementioned enhanced smoothing properties play a major role in establishing 
the well-posedness of the GCSIE in the smooth case. We address in this paper the issues of well- 
posedness and well-conditioning of the GCSIE formulations in the Lipschitz case by making use 
of deep results from harmonic analysis m and certain duality pairings. We also present a direct 
counterpart to the GCSIE formulations which we refer to as Regularized Combined Eield Integral 
Equations (CFIER). From the numerical point of view, the main advantage of the direct CFIER 
formulations is given by the fact that they employ as unknowns Dirichlet and Neumann traces of 
transmission problems, whose singular behavior around corner points is well understood m and 
thus can be resolved by graded meshes towards corners. This brings us to the second major point 
of this paper, high-order Nystrom discretizations of transmission boundary integral equations in 
two-dimensional Lipschitz domains. 

High-order Nystrom methods typically employ graded meshes in order to deal with singulari¬ 
ties associated with solutions of boundary integral equations in domains with corners |22l[I7l[l6|. 
One issue that arises in this regard is the possibly unbounded nature of such solutions in the 
neighborhood of corners in the case of integral formulations of the second kind. While in several in¬ 
stances the issue can be avoided by resorting to alternative integral equation formulations [2| whose 
solutions are regular enough (e.g. Holder continuous), in many others, including the case of trans¬ 
mission problems in domains with corners, the unboundedness of solutions cannot be avoided. Two 
main approaches to tackle the unbounded nature of solutions of integral equations of the second 
kind have been recently introduced: (a) one that relies on incorporation of the known asymptotic 
infinite behavior of solutions in the vicinity of corners and exact cancellations of infinite quanti¬ 
ties |9|; and (b) one that uses jacobians associated with graded meshes in order to introduce more 
regular weighted solutions as new unknowns of newly weighted integral formulations of the second 
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kind (HU El [ISl HH EH] . We pursue a version of the latter approach in this paper in conjunction with 
simple modifications of a Nystrom method based on global trigonometric interpolation, singular 
kernel-splitting, and analytic evaluations of integrals that involve products of certain singular func¬ 
tions and Fourier harmonics |26l E9] . Several ideas in this approach were introduced in |28j for the 
discretization of the first kind boundary integral equation formulation of transmission problems |13j 
in domains with corners. The method incorporates sigmoid transforms |22j within parametrizations 
of domains with corners and it uses the Jacobians of these transformations as multiplicative weights 
to define new unknowns. Specifically, the focus of the paper is on direct integral formulations of 
transmission problems whereby the unknowns are the Dirichlet and Neumann traces of solutions 
of transmission problems on the Lipschitz boundary. A weighted Neumann trace defined as the 
product of the derivatives of the sigmoid parametrizations and the usual Neumann trace of solu¬ 
tion of transmission problems is introduced as a new unknown; given that the derivatives of the 
parametrizations that incorporate sigmoid transforms vanish polynomially at corners, the weighted 
traces are more regular for large enough values of the order of the polynomial in the sigmoid 
transform. Introducing new weighted unknowns also require definition of new weighted boundary 
integral equations that involve weighted versions of the four scattering boundary integral operators. 
It turns out that the kernel splitting techniques originally developed for smooth curves |23j can 
be easily extended to the weighted boundary integral operators, delivering a high-order Nystrom 
discretizations for the various formulations considered in this paper. We give ample numerical evi¬ 
dence that in the high-contrast, high-frequency regime the single integral equation formulations and 
our novel regularized formulations have superior spectral properties over the classical formulations 
of transmission problems, giving rise to important computational savings. 

The paper is organized as follows: in Section]^ we formulate the Helmholtz transmission prob¬ 
lem we are interested in; in Section we recount the definition of the four scattering boundary 
integral operators and we discuss several boundary integral formulations of the Helmholtz trans¬ 
mission problem; in Section we establish the well-posedness of several of the boundary integral 
formulations discussed in this paper, including the regularized formulations CFIER; finally, in Sec¬ 
tion we present high-order Nystrom discretizations of the various boundary integral equations 
considered in this paper and we carry out a comparison between solvers based on this formulations 
that emphasizes the benefits that can be garnered from the use of single integral equations and the 
regularized integral equations. 


2 Integral Equations of Helmholtz transmission problems 


We consider the problem of evaluating the time-harmonic fields and that result as an incident 
field impinges upon the boundary F of a homogeneous dielectric scatterer D 2 which occupies a 
bounded region in We assume that both media occupying D 2 and its exterior are nonmagnetic, 
and the electric permitivity of the dielectric material inside the domain D 2 is denoted by 62 while 
that of the medium occupying the exterior of D 2 is denoted by ei. The frequency domain dielectric 
transmission problem is formulated in terms of finding fields v} and v? that are solutions to the 
Helmholtz equations 


= 0, in D 2 , 

Au^ + kfu^ = 0, in Di = \ 1 ) 2 , 
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where the wavenumbers ki,i = 1,2 are defined as ki = ijjy/el, z = 1 , 2 in terms of the frequency oj. 
The incident field is assumed to satisfy 

An®®®" + = 0 in^a, (2.2) 

where D 2 is an open neighborhood of Ha- Therefore n®”" has to be smooth, actually analytic, in T^a, 
which includes plane waves or spherical waves from a point source placed in the exterior domain 
Di. 

In addition, the fields v}, n®®®", and v? are related on the boundary T by the the following 
boundary conditions 


+ y^n®®®" = yfjn^ on T 

yj^n^ + 7 Arn®®®" = py^n^ on T (2.3) 


with p > 0. In equations (2.3) and what follows 7 ) 5 , z = 1,2 denote exterior (z = 1) and interior 
Dirichlet traces (z = 2). Similarly 7)^,z = 1,2 denote exterior and interior Neumann traces taken 
with respect to the exterior unit normal on T. When both Dirichlet, respectively Neumann, traces 
coincide, we will simply write 7 ^), respectively 7 Ar. (Notice that 7 ^zz®®®" = y^jzz®®®", 7 )(fZZ®®®" = y^zz®®®" 
and therefore we are allowed to use in ( |2.3[ )). 

We assume in what follows that the boundary T is a closed Lipschitz curve in Depending 
on the type of scattering problem, the transmission coefficient p in equations (2.3) can be either 
1 (E-polarized) or ei/e 2 (H-polarized). We furthermore require that v} satisfies the Sommerfeld 
radiation conditions at infinity: 


lim {dv} / dr — ikiv^) = 0 . 

|r|^oo 


(2.4) 


(Here d/dr is the radial derivative). Note that under these assumptions the wavenumbers ki, 
z = 1,2 are real numbers. It is well known that in this case the systems of partial differential 
equations (2.1)-(2.2) together with the boundary conditions (2.3) and the radiation condition (2.4) 


has a unique solution [25l [21]. The results in this text can be extended to the case of complex 
wavenumbers ki,i = 1 , 2 , provided we assume uniqueness of the transmission problem and its 
adjoint, that is, for the same transmission problem but with wavenumbers ki for D 2 and k 2 for Di 


3 Boundary integral formulation for the transmission problems 


A variety of well-posed integral equations for the transmission problem (2.1)-(2.3) exist [siKiaizii 
0]. On one hand, integral equations formulations for transmission problems can be formulated as a 
2 x 2 system of integral equations which can be derived from (a) Green’s formulas in both domains Di 
and D 2 , in which case they are referred to as direct integral equation formulations |13[l21j . (b) from 
representations of the fields zz-^, j = 1 , 2 in forms of suitable combinations of single and double layer 
potentials in both domains Di and D 2 , in which case they are referred to as indirect integral equation 
formulations [25], or (c) from Green’s formulas and suitable approximations to exterior and interior 
Dirichlet-to-Neumann operators, in which case they are referred to as regularized combined field 
integral equations or generalized combined source integral equations |3]. On the other hand, integral 
equations formulations for transmission problems can be formulated as single integral equations 
which can be derived from (d) Green’s formulas in one of the domains and (indirect) combined field 
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representations in the other domain [2T]. The strategies recounted above lead to Fredholm second 
kind boundary integral equations for the solution of transmission problems [251 [mill [5], at least in 
the case when the curve F is smooth enough {C^ suffices). The first part of this paper is devoted 
to establishing the well-posedness of the boundary integral equations of the type (c) and (d) in the 
case when the curve F is Lipschitz. To this end, we begin by reviewing the definition and mapping 
properties of the various scattering boundary integral operators. 

3.1 Layer potentials and operators 

We start with the definition of the single and double layer potentials. Given a wavenumber k such 
that ?R.k > 0 and ^k > 0, and a density (p defined on F, we define the single layer potential as 

[SLk{(p)]{z) := Gk{z - y)(p{y)ds{y), z e \ F 

and the double layer potential as 

[DLk{^p)]{z) := z e \ r 

where Gfc(x) = jFrQ^^(/c|x|) represents the two-dimensional Green’s function of the Helmholtz 
equation with wavenumber k. The Dirichlet and Neumann exterior and interior traces on F of the 
single and double layer potentials corresponding to the wavenumber k and a density (p are given by 

loSLkip}) = -f]:,SLk{ip) = -fnSLkiip) = Skip 

-flfSLkiip) = + kJ<p j = l,2 

^ (3.1) 

y^^DLk{ip) = {,-\y^^\^Kkip j = l,2 

lljDLkiip) = lljDLkiip) = -fNDLk{ip) = ip. 

the operators Kk and , usually referred to as double and adjoint double layer 
operators, are defined for a given wavenumber k and density ip as 

{Kkip)ix) := ip{y)ds{y), x on F (3.2) 

and 

{Klip){x) := <fiy)ds{y), x on F. (3.3) 

Furthermore, for a given wavenumber k and density ip, the operator Nk denotes the Neumann trace 
of the double layer potential on F given in terms of a Hadamard Finite Part (FP) integral which 
can be re-expressed in terms of a Gauchy Principal Value (PV) integral that involves the tangential 
derivative dg on the curve F 

(JV.^Xx) := FF 

= Gk{x - y)(n(x) • n(y))(/?(y)ds(y) + PV ^ dsGk{x - y)dsip{y)ds{y). (3.4) 


In equations (3.1 
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Finally, the single layer operator Sk is defined for a wavenumber k as 


{Skip)ix) := j Gk{x-y)(p{y)ds{y), xonT (3.5) 

for a density function tp defined on F. 

We commit here another slight abuse of notation and denote in what follows SLi,DLi, Si, Ni..., 
for layer potentials and operators corresponding to ki. That is, 


SLi = SLk, DLi = DLk„ Si = Sk„ etc.. 

Analogously 

SL 2 = SLk, DL 2 = DLk„ S 2 = Sk„ etc.. 

We stress that the context will avoid any possible confusion between indices and wavenumbers. 
This convention helps us to enhance and lighten the notation. Notice that Green identities can be 
now written in the simple form: 


= (-lySLji'ylrU^) - {-1)W 


(3.6) 


Similarly, 


G- — f 
~ 2 


I 


+ (-iy 


-Kj S, 
-N, K], 


j = l,2 


(3.7) 


are the exterior/interior Calderon projections associated to the exterior/interior Helmholtz equa¬ 
tion: 








(3.8) 


We recall that from (3.7)-(3.8) one deduces easily 


SkNk = -\l + Kl, NkSk = -\l + {Klf, NkKk = KjNk, KkSk = SkKj. (3.9) 


We end this section noticing that since solves Av + kfv = 0 in D 2 , it holds 


and therefore 


Cl 


IDU 

INU'-' 


= 0 . 


- SLi (77V«'”") + DLi = 0, in 


(3.10) 

(3.11) 


3.2 Boundary integral equations 

Let us introduce first the total field given by 


u* = 


V? 


in Di 
in D 2 


The unknowns in the direct formulations we consider in this paper are 
IDU^ = iW = 7Ar(«^ + 


(3.12) 


(3.13) 
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the Dirichet and Neumann data for the total held from the unbounded domain. We could then 
construct using Green’s identities (|3.6|). Actually, for the exterior solution it holds as well 


u 


^ = -SLii'ylfU^) + DLi{jdu^) 


(3.14) 


which is consequence of the dehnition of the total wave n* and of the identity (3.11) which ensures 
that the contribution from u*"'’ in the representation formula cancels out. 

We hrst have the following hrst kind integral equations due to m with a positive dehnite 
principal part 


CFK 




iW. 



-(i^i+ATa) 

-{N^+pN2) 


-1 


(P 

{Kj + KJ 


S2 + Si) 

2 ) 



IDU^ 


ryj^yinC 


fiW. 




(3.15) 


This system of integral equations can be easily derived from (3.7)-(3.10) and the boundary condl 


tions (2.3). Obviously, these integral equations are understood in the a.e. sense. We refer to the 


system of boundary integral equations (3.15) as CFIEFK for combined field integral equation of 
first kind. 


A widely used hrst kind boundary integral formulation of the transmission problem (2.1)-(2.3) 
is due to Kress and Roach cf. 


and it consists of the following pair of integral equations: 


CSK 


iDvf 

/ P“' + 1 

7 

+ 

1nu\ 

■■ 2 

I 


(K2 - p-^Ki) 
-{N 1 -N 2 ) 


P-HS1-S2) 

(Kj-p-^Kj) 


1 ) 

IDU^ 



\ 

1nu\ 




(3.16) 

In what follows we refer to the integral equations (3.16) as CFIESK (combined field integral equation 
of second kind). 

We introduced recently Generalized Combined Source Integral Equation (GCSIE) formulations 
of transmission problems HE] and we established their well-posedness in the case when F is regular 
enough. We consider here a direct counterpart of the GCSIE formulations which can be obtained 
by combining the previous formulations in the form 


p+l 


CSK + 


l + p 




-pN^ 


CFK j 


7£)U* 

1 

I 

25,1 

ryj^yinC 


P + l 

[-2piV, 

Pl\ 



(3.17a) 


Here, k is a complex wave number with positive imaginary part. This parameter can be appropri¬ 
ately taken to improve the spectrum of the underlying operator. This leads to faster convergence, 
when discretized, of Kyrlov iterative methods, such as GMRES, for the linear system (see subsection 
5.4). Dehning 


R := 


1 


P + l 


I 

-2pN^ 


25, 

pi 


we observe that (3.17a) can be written as 


GFK 




ri/ + K 2 

pN2 


-p-^S2 

ll-Kj 


iDvf 

iW. 


+ RCFK 


IDU" 

llu\ 


= R 


iNVd' 


(3.17b) 


We will refer to this formulation as Regularized Combined Field Integral Equations (CFIER). 


Another possible formulation of the transmission problem (2.1 )-(2.3) takes on the form of single 
integral equations [2T]. The main idea is to look for the field as a single layer potential, that is 


n^(z) = -2[5L2 P](z), z G D 2 , 
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where /U is an unphysical density defined on F. The transmission boundary conditions (2.3) and the 
trace relations (3.1) imply first 


IDU^ = -252^, = -p{I + 2 K 2 )p 




and second, from (3.14), 


tt^(z) = pSLi[{I + 2iF2^);u](z) - 2DLi[S2p]{z), z E Di. 

Using these representations of the fields and the Neumann and Dirichlet traces of and on 
r are used in a Burton-Miller type combination of the form = 

—cf [To] to lead to the following boundary integral equation 


where 

and 


SIE p := -~ ipSp = 0 7 ^ r/ E 

K = -K^{pi - 2Kj) - pKj {I + 2Kj ) + 2(iVi - N2)S2 
S = -pSi{I + 2Kj ) - (/ - 2iFi)52, 


(3.18) 


We refer in what follows to equation (3.18) as SCFIE. The coupling parameter p in equa¬ 


tions (3.18) is typically taken to be equal to ki. 


4 Existence and uniqueness for the boundary integral formula¬ 


tions CFIEFK (3.15), CFIESK (3.16), CFIER (3.17), and SC¬ 


FIE (3.18) 


In this section we state the well-posedness of the boundary integral formulations on Lipschitz curves 
presented in previous section. Let us point out that, with very minor and direct modifications, the 
proofs of the main results of this section can be adapted to the case of Lipschitz domains in 
3D. Hence, these formulations are well-posed as well, and the associated integral operators satisfy 
the same properties in the Sobolev frame. Since in this paper we have focused our attention to 
bidimensional domains, we have preferred to work only on Lipschitz curves for the sake of simplicity. 

One of the basic tools we will use in the analysis developed in this section consists in comparing 
the boundary operators for the Helmholtz equation with those for Laplace operator, Sq, Kq, Kq , 
and Nq which are defined as the corresponding boundary Helmholtz operators in section 3 using 

Go(x) = -^log|x|, 

the fundamental solution of the Laplace operator, instead. 

4.1 Sobolev spaces and functional properties of boundary operators 

For any H C domain with bounded Lipschitz boundary T, we denote by H^{D) the classical 
Sobolev space of order s on H (see for example in the excellent textbooks [301 Ch. 3] or [H 
Ch. 2]). We consider in addition the Sobolev spaces defined on the boundary T, H^{T), which 












are well defined for any s G [—1,1]. It is well known that for smoother T the range for these s 
can be widened but we do not make use of these spaces in this section. We recall that for any 
s>t, C S G {Di, D 2 ,T} with continuous and compact embedding. Moreover, and 

(l/*(r)) = when the inner product of H^{T) = L‘^{T) is used as duality product. 

It is well known that 'yj^ : — )• H^{T) is continuous for s G (0,1) and if 

HkiDj) ■■= {U G H^{Dj) : AU G L^Dj)} , 

endowed with its natural norm, then 77 V : H^{Dj) —)• is continuous for s G (1/2,3/2). 

The space H^(r), and its dual are then the limit case from several different perspectives. 

Let Vr be the tangential derivative defined as 

Vr(t/|r) = 5.(17|r)r 

where r denotes (one of) the unit tangent vector fields to T, and dg the tangential derivative with 
respect to r. It is known that an integral expression for the (or an equivalent) inner product in 
H^(T) is given by 

{u,v)ffi(^py.= j Vrrt • Vru + J uv = (Am, An), AM:=VrM + Mn (4.1) 

We commit here a slight abuse of notation, and denote with the same symbol (•, •) the non-complex 
integral inner product in (L^(r)) . 

Theorem 4.1 Let D 2 be a hounded domain, with Lipschitz boundary T, and set Di := \ D 2 . 

Then, if X ^ C°°(M^) with compact support, 

. SLk : ^ H^+^/^D2), xSL^ : F^(r) ^ 

• DLk : iL"+3/2(D2), xDLk : 

are continuous for s G [—1,0]. Moreover 

• Sk : H%r) iL"+i(r) 

• Kk : ^ F"+i(r) 

. Kj : H^{T) ^ H^{T) 

• Nk : iL'+^T) ^ H%T) 
are continuous for s G [—1,0]. 

A proof for the intermediate values s G (—1,0) can be found in [l2l Th 1] (see also [30l Th 6.12]). 
The proof for A; = 0, the Laplace operator, and s = —1,0 can be found in m (see also the comments 
following Th 6.12 in [30j 1. For k ^ 0 the argument follows by showing that the difference between 
the corresponding Laplace and Helmholtz boundary operators are smooth enough. We refer to the 
discussion at the end of Chapter 6 in [30] or that following Theorem 1 in [T2| and references therein. 
For the sake of completeness, we will give next a proof of this result. This result will be actually 
used to prove the well posedness of the different formulations considered in this paper. 

We will need first this technical lemma which will be proven for the sake of completeness. 
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Lemma 4.2 The integral operators 

A(-K'a:i - Kk^), Nk^ - Nk2 

have weakly (integrable) singular kernels. 


Proof. Without loss of generality we can assume in this proof that /c 2 = 0. The kernel of 
A{Kki — Kq) is then given (a.e.) by 

-n(x) [(V(Gfci - Go)(x - y)) • n(y)] - t(x) t’^(x) - Go)(x - y)) n(y) , x,yer. 

( above denotes the Hessian matrix of G), whereas 

-n'^(x) - Go)(x - y)) n(y) 

is the kernel of — Nq. 

Then it suffices to show that for any R > 0 there exists Cr > 0 such that 

|V(Gfc, -Go)(x)| + ||v2(Gfc^ -Go)(x)|| <GR(l + log|x|), for 0 < |x| < i?. (4.2) 


where || • || above is any matrix norm. 

The proof of this bound relies on analytical properties of the Hankel functions, namely, the 
behavior at zero and appropriate decompositions of these functions. We point out that a deeper 
analysis on this topic will be carried out in Section so that we limit ourselves the exposition of 
the properties we are going to use. 

With the identity we can obtain first 

V(Gfci - Go)(x) = (^-^H[^\ki\:>c\)ki\Dc\ + (4.3) 

(we follow the convention of writing the gradient as a row vector). Note now that 

^zH[^\z) = -^Ji{z)z\ogz + ^ +di{z)z‘^ 

where Ji is the Bessel function which is known to be smooth, with Ji(0) = 0, J((0) = 1/2 and di 
being smooth as well. From this decomposition we deduce easily that V(Gfc^ — Go) satisfies 
On the other hand, using that 

[zH[^\z)y = zHI,^\z). 

we can easily show that 

V2(Gfci - Go)(x) = -*^iLj^^(A;i|x|)|^xx'^ + ^JiLi^^(A:i|x|)A:i|x| - ^|^xx^- 

where / is 2 x 2 identity matrix. 

Since in addition 

= -^Jo{z)logz + co{z) 

where Jq (the Bessel function) and cq are smooth, a simple inspection shows that the entries of 
the Hessian matrix can be bounded by Glog(|x|) on any punctured ball around 0 with appropriate 
constant G (which depends obviously on ki and on the diameter of the domain). ■ 


(4.2). 
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Theorem 4.3 Let ki ^ k 2 - Then 
. Sk,-Sk,:H-\T)^H^{T) 
. Kk, - Kk, ■. H^T) ^ H\T) 

. Nk,-Nk,:H\T)^H^{T). 
are continuous and compact. 


Proof. It is a well known consequence of Lax Theorem (see |24[ Th 4.12]) that integral operators 
with weakly singular kernels define compact operators in L^. From this fact, and Lemma 
see that 

Nk, - Nk, : H\T) ^ H^{T) 

is compact. 

Similarly, 

A{Kk,-Kk,):H\T)^[H^{r)]^ 

is also compact (again from Lemma |4.2| ). 

Consider now a weakly convergent sequence {un)n- Then (^A{Kk^—Kk 2 )un)^ converges strongly, 

by the compactness of the operator, in [iL°(r)]^ and, because ( |4.1| ), [{Kk^ — Kk 2 )un)^ is strongly 
convergent in H^{T). In other words, Kk^ — Kk 2 '■ H^{T) H^iT) is compact as well. 

A duality argument proves now the result for . 

Finally, assume without loss of generality, that ki is taken so that : H^{T) —)■ H~^(T) is 
invertible. (It suffices to take as ki a pure imaginary number). Then 


4.2 


we 


Sk,-St, = N;^"lNi,^St,-Nt,St,]+N;:^^(Nk,-Ni„)Si„ 

= - Tlf] + 

= "if Tl - KlTl + [< - (4-4) 


where we have applied the second identity in (3.9). In view of 
the operators involved, we conclude that Sk 2 — *S'fci : H~^{T) - 
The proof is now finished. ■ 


[ |4.4| ) and the mapping properties of 
■ H^iT) is continuous and compact. 


The last ingredient in our proof is this result due to Escauriaza, Fabes and Verchota [13]. In 
this result, Kq, Kq are the double and adjoint double layer operator for Laplace equation (which 
obviously correspond to A: = 0). 


Theorem 4.4 For any Lipschitz curve and A 0 (—1/2,1/2], the mappings 

XI+ Ko: H%r) LC(r), \I + K^ : H-%T) H-^T) 

are invertible for s G [0,1]. Moreover 

\l + Ko: H^{r) ^ h + ^ H-^{T) 

are Fredholm of index 0. 
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Proof. In |37t Th 3.1 and 3.3] it is proven that —\l + Kq : H^(T) —)■ H^{T) is invertible 
for s = 0,1. (Note that in this paper ^ log |x| = —Go(x) is taken as fundamental solution.). By 
interpolation of Sobolev spaces and a transposition argument, we can prove the result for A = —1/2. 

In [m Th 2] the result is proven for |A| > 1/2 and s = 0. To extend the result for the remaining 
s we can adapt the argument outlined in I37j. Suppose that Sq : H^{T) —)• H^{T) is invertible. 
Then, So{XI + Kq)Sq^ : H^{T) —)• H^{T) is invertible as well. But Calderon identities ( |3.9[ ) imply 

SoiXI + Ko)Sq^ = {XI + k;[ )SoSo^ = (XI + ) (4.5) 

from where we conclude the invertibility, for A 0 (—1/2,1/2], of XI + KJ : 77® (T) —>■ 77®(T) first 
for s = —1 and next, by interpolation, for s G [—1,0]. Again, by transposition, we can show that 
XI + Kq : 77® (r) —7- 77® (r) is invertible for s G [0,1]. This argument breaks down when So is not 
invertible, that is, when the logarithmic capacity of the curve T is 1. For this case we consider 

So<f := So(p + y 7’- 

It can be shown that .So : 77®(r) —k 77®+^(r) is continuous and invertible for s G [—1,0] (see [HUl 
Ch. 8]). Moreover, since 

Ik„, = ks{1,) 

we still have KqSq = SoUg, which makes possible to extend the argument in 
to this case as well. 

Finally, in [371 Th 4.2] it is proven that ^I+KJ : TQ(r) —>■ Tq(F) is invertible where 7yQ(F) is the 
subspace of functions in T^(F) = 77’^(r) which zero mean. Then dim A^(^7 + 77 q~) = codim R(\l + 
Kq ) = 1, which in particular implies that ^I + Kq \ H^{T) —)• 77*^(r) is Fredholm of index zero 
and, by transposition, so is ^7 + Kq \ 77^(F) —?■ 77‘^(r). The same argument as before extends this 
result to 77®(r). ■ 

Observe that in the proof of this last result we have assumed that F is simply connected. 
Otherwise dimA^(|7 + 77) = codim 77(^7 + 77 q~) equals to the number of simply connected 
components, and the argument is still valid with this very minor modification. 


(4.5), with So instead. 


4.2 Well-posednees of the boundary integral equation formulations CFIEFK (3.15), 


CFIER (3.17), and SCFIE (3.18) 


In this section we state and prove the well-posedness of the integral equations formulations for the 


transmission problems. For Costabel-Stephan CFIEFK (3.15) and Kress-Roach CFIESK (3.16) for¬ 
mulations, the stability has has been already proved, but only in the space 77^/^ (F) x 77-1/^) da 
36] - We will however extend these results to a wider range of Sobolev spaces. To this end, we make 
use of the skew-symmetric bilinear form 


((/, 7), (ff, ^)) ■■= J^f^- 57, (/, 7), (5, i’) G R'/'(r) X 77-1/2 (D 

which gives a non-usual representation, but more convenient for our purposes, of the duality product 
between 77i/^(F) X 77 i/^(F) and itself. Obviously, the integrals above have to be understood in a 
weak sense if (p and ■0 are not sufficiently smooth. 
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We point out that if A : x H ^ x H with A = then 

the transpose operator 


{A\f,ip),{g,^)) := ((/,(^),A(5,V’)), V(/, (<7, G ^'/'(F) x 


is given by 


A* = 


1' 

A^ 

1' 


-1 

lb. 

1- 

1 

-1 

-1 


1- 

1 

1 


(4.6) 


Here, Ah is simply the adjoint of Aij in the bilinear form defined by the integral product in F. 
Observe that the familiar identity (AH)* = H*A* still holds. 

Theorem 4.5 The following operators 


CSK : H"(F) X iH-*(F) ^ H%T) x H'*-*(F), CFK : H"(F) x H"(F) x H^-^T) 

are invertible with continuous inverses for all s G [0,1] 


Proof. Assume p 7 ^ 1- Note that from (3.16) 

p-1 


CSK = 


P 


p+i 

2(p-l) 


I + Ko 


p+i 

2(p-l) 


I+ K, 


T 


CSKo 


+ 


K 2 - Ko - p-^Ki - Ko) p-^(Si-S2) 

-(W - N 2 ) Kj - Kj - p-^iKj - Kj ) 


(Kq, Kq denote here the double layer and adjoint double layer operator for the Laplace equation). 
Clearly, 

[- 1 / 2 , 1 / 21 . 


By Theorem 4.3, Li : H^{T) x H~^(T) —)■ H^{T) x H^{T) is continuous and compact. On the other 
hand, CSKq : H^{T) x —)• H^{T) x H^“*(F) is invertible by Theorem 4.4, Therefore, CSK 

is Fredholm of index zero with kernel in H*(F) x H**(F). For p = 1 the proof is even simpler since 
Theorem 4.3 shows now that CSK is a compact perturbation of the identity in H^fT) x H"-*(F) 
and the same argument can be applied to prove that the kernel is contained in H*(F) x H^(T) as 
well. 

The injectivity can be shown using Green identities and the uniqueness of the transmission 
problem cf. m- We emphasize that the regularity of the kernel of CSK is crucial to make the 
arguments in m valid. 

For the second operator, we first note that 


CFK = 


-2Ko 

-(l + p)A^o 27^0^] 

CFKo 


-I-L2 
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where L2 : H^{T) x H ^(F) —)■ H^{r) x H^(T) is again compact. Calderon identities in (3.9) yield 
that (see also (4.6)) 


CFK(,CFKo = 


2Ko 


-‘^S„ 


-2K„ 


‘fSo 


_(1 + p)No -2Kj \ [-(1 + p)No 2Kj 

-4K^ + ^ SoNo ^ (KoSo - SqKJ ) 
-2{p+l){NoKo-KjNo) -4{Kjf + 


{p+^y 


NoSo 


(p-lf 


(p+1)" 

4(p-l)2 


I + K^ 


(p+i)^ j \ (T-r 


4{p-l) 




jp-l) 

P 


2 r 


{-f3I + Ko){/3I + Ko) 


i-/3I + Klm + K^)^ 


/3:= 


P+1 

2p-2' 


Theorem 4.4, observe that |/3| > 1/2, implies that this mapping is invertible, and therefore, so is 
CFKq. We then conclude that CFK : H^{T) x H^{T) x is Fredholm of index 

zero . If p = 1, 


CFK^ CFKn = - 


I 


I 


which makes the argument still valid. The injectivity, which implies the invertibility by Fredholm 
alternative, follows using classical arguments cf. [36]. ■ 


We would like to point out that the original proof for the invertibility of CFK is based, in a 
nutshell, on showing that its principal part is coercive (strongly elliptic) in the space x 

^-i/2(r) Section 5]. 


Theorem 4.6 Assume that the wavenumber k in (3.17) has a positive imaginary part. Then the 
operators 

CFK : FC(r) X H%T) x 

are invertible with eontinuous inverses in the spaces H^{T) x 77*“^(F) for all s G [0,1]. 

Proof. We first prove that GFK is a compact perturbation of an invertible operator. Starting 
from (3.17b), and proceeding as above, we derive 
r , .-ic_ 1 1 


GFK = 


1/ + K0 -p-^So 

pNo 


1/ _ /VT 
2-' -^0 


+ 


P+1 


I 2So 

-2pNo pi 


CFK0 + L3 =: GFK0 + L3. 


Lq Ro 

with L3 : H^{T) x H~^{T) —)■ H^{T) x H^{T) compact. Straightforward calculations show that 


LqLq — 0 


CFK() R(, Ro CFKo = 


4(P-1)^ 
(P + 1)2 

:fKo4 

2(P-1) 

P+1 


■^0(-2^^ + ^0)(2^/ + i^0) 




L*o Ro CFKo + CFK*o R(, Lo = 

■(-/ + 2K2)(^/ + Ko) 


{-I + 2{Kjr) {^^I + iKj)) 
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Therefore, 


where 


Note that 


GFK*o GFKo = 


Ai 


^7 


(4.7) 


Ai : = 


4(P-1)- 

{p+l? 


Ko + 


P + ^ ('^3 _ 3(/3+ 1) j^2 , /5 + 1 


2{p-l) 


K^ + 


2{p-l) ^ 2{p-l) 


= 77^ ( 2 ^^^ + ^ 0 ) ( - Ai(p)/ + Ko){- \2{p)I + K,){- X,{p)I + Ko) 
where Xj{p) are the roots of 

= ^3 _ + P +1 


2{p-l) 2{p-l) 


which can be shown to be all real, and lying outside of [—1/2,1/2] for any p G (0, 00 ) \ {1} cf. 


Lemma 4.7 below. Theorem 4.4 implies the invertibility of Ai and therefore of GFKq. 


For p = 1 the proof clearly breaks down, but it can be studied separately and deduce that (4.7) 
still holds with 


1 


1 


which again is invertible. 

In short, we have shown that GFK : H^{T) x —)• F7^(F) x for s G [0,1], is 

Fredholm of index zero. To prove it is invertible we will consider the adjoint operator given by (see 


(3.17b)) 


GFK* = 


+7-772 

1 

7 

+ 

'{Ki + K 2 ) 

-{p-^S2 + Si) 

1 

■ pi 

-2S^ 

-pN2 

ll + Kj 

Ni + pN2 

-{Kj + k:i)_ 

1 + p 

2pN^ 

I 


We note that the operator GFK* is exactly the operator corresponding to the GGSIE formula¬ 
tions HE]- Glearly 

GFK* = GFK[, L*3 

where GFKf, : H^{T) x -G H^{T) x is invertible for s G [0,1] and L* : H°(T) x 

i7-*(F) ^ H^{T) X 77**(F) is compact. Then, the null space of this operator is contained in 
77* (F) X 77** (F). We can then follow the arguments in [5] and conclude GFK* is injective (the 
arguments are still valid for Lipschitz domains once we ensure that the elements of the null space 
are smooth enough), and therefore, from Fredholm alternative, invertible from where one derives 
the invertibility of our operator GFK. ■ 


Lemma 4.7 For any positive p 7 ^ 1, the roots of the polynomial 

3 _ 3 (p + 1) 2 , P + 1 
2 (p-l) ^ 2 (p-l) 


are all real and lie outside the interval [— 1 / 2 , 1 / 2 ]. 
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Proof. We can reduce the problem to study the polynomial 


q/ 3 {x) = — 3/?x^ + /3 

ior 13 = ^ GR\[-1/2,1/2]. 

Assume /3 > 1/2. Then (^^(±1/2) = ±| + j. That is, ^^^(±1/2) > 0. Also, g/3(2/3) = 
4/3 — /3^) < 0, and thus qj 3 has a zero between 1/2 and 2/3. Furthermore, since lima;^ipoo qjsix) = 

=Foo, it follows that qp must have one (real) zero in (—oo,—1/2), and another zero in (2/3, oo). 
Hence, we have shown that the three roots of the polynomial are real and he outside the interval 
[-1/2,1/2]. 

For /3 < —1/2, one can show, similarly, that has three real roots, two of them smaller than 
—1/2, and the other one larger than 1/2. ■ 


We end this section studying the SCFIE formulation (3.18) which again turns out to be stable, 
although in weaker norms. 


Theorem 4.8 The operator SIE : H^{T) 
invertible for any s G [—1,0] 


H^{T) associated to the formulation SCFIE (3.18) is 


Proof. It can be easily seen that if p 7^ 1 we have 


SIE = 2(1 -p)( - J- 


2 { 1 + (- 


4(1-p) 1-p 

p+i 
'2(p-l) 


^0 + [K^ 


+ L4 


I + Kq ) + L 4 


where L4 : —)• can be checked to be compact. We have that the operator + Kq is 

Fredholm of index zero (cf. Theorem 4.4) while the operator — 2(^^\) ^ invertible, and 

thus the main part of the SIE operator is kredholm of index zero as well which is enough for our 
purposes. Fredholm alternative implies now the invertibility of SIE. (We refer to the seminal paper 
for a proof of the injectivity). The case p = 1 can be treated similarly. ■ 


5 


High-order Nystrom methods for the discretization of the 
mulations CFIEFK (|3l^, CFIESK (|3l^, CFIER (|M7|), 
SCFIE (|3l^ 


for- 

and 


We present in this section Nystrom discretizations of the formulations CEIEEK (3.15), CEIESK (3.16), 
CFIER (3.17b), and SCFIE (3.18) presented in the previous section. The discretizations of some of 
these formulations, for smooth curves, have been already analysed in [7]. The key component for 
polygonal domains is to use sigmoidal-graded meshes that accumulate points polynomially at cor¬ 
ners and to reformulate the aforementioned systems of integral equations in terms of more regular 
densities and weighted versions of the boundary integral operators of Helmholtz equations. 
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5.1 Weighted boundary integral operators for Helmholtz equation 


We assume that the domain D has comers at xi,X2,... ,xp whose apertures measured inside D 
are respectively 61 , 62 ,... ,9p, and that F \ {xi, X2 ,xp} is piecewise analytic. Let {xi{t),X 2 {t)) 
be a 2 tt periodic parametrization of F so that each of the curved segments [xj,Xj_|_i] is mapped 
by {xi{t), X 2 {t)) with t G [Tj,Tj^i]. We assume that xi{t),X 2 {t) are continuous and that on each 
interval [Tj,Tj^i] are smooth with + {x 2 {t))‘^ > 0 (the one-sided derivatives are taken for 

t = Tj,Tj^i). Consider the sigmoid transform introduced by Kress in |22j 


w{s) 

v{s) 


Tj+i[v{s)]P + T,[l-vis)]P 
[u(s)]P + [1 -u(s)]P 

(p - 2 ) 


Tj < s < Tj+i, I <j < P 

1 2s — Tj — Tj^i 1 
p Tj+i — Tj 2 


(5.1) 


where p > 2. The function re is a smooth, increasing, bijection on each of the intervals [Tj,Tj^i] 
for I < j < P, with w^^\Tj) = w^^\Tjj,.i) = 0 for 1 < /c < p — 1. We then define the new 
parametrization 


X{t) = {xi{w{t)),X2{w{t))) 


extended by 27r—periodicity, if needed, to any t G M. 


A central issue that collocation discretizations of the integral equations CFIEFK (3.15), CFIESK (3.16), 


and SCFIE (3.18) is confronted with is the possibly unbounded nature of their solutions at corners. 


In what follows we deal with this issue by introducing new weighted unknown densities defined as 

:= (7iv^‘)(x(t)) |x'(t)| (5.2) 


for equations CFIEEK (3.15) and CFIESK (3.16), as well as a new weighted unknown density 
dehned as 

^"'(t) ;=//(x(t))|x'(t)| (5.3) 


for the equation SCFIE (3.18). 


According to the classical theory of singularities of solutions of elliptic problems in non-smooth 


domains, the solutions of the integral equations CFIEFK (3.15), CFIESK (3.16), CFIER (3.17), 
and SCFIE ( |3.18 ) exhibit corner singularities |13j . In the case of smooth incident helds (e.g. 
pla ne w ave incidence) we have that G 77^(F) x L^(F). Given that by Theo¬ 

rem 4.5 the operators CSK (and/or CFK) are invertible in the space F7^(F) x L^(F) we obtain that 
G H\T) X Lp'iT). Consequently, it follows from Sobolev embedding theorems that 


7 p)?x* is Holder continuous on F. Also, given that by Theorem 


4.8 


the operators SIE are invertible in 


L^(F), we obtain that the solution p of the SCFIE equation (3.18) belongs to L^(F). In conclusion, 
in the case of smooth incident fields the integrable corner singularities of yAru* and p are 
mollihed by the corner polynomially vanishing weights |x'| for large enough values of the exponent 
p of the sigmoid transform. Indeed, for large enough values of p the 27r periodic functions 
and p'^ vanish at Tj for all 1 < j < P and are regular enough. 


Remark 5.1 More precise statements can be made about the nature of the functions yjy’^u* and 
p^ if we resort to the corner asymptotic behavior of Indeed, it was shown in that 
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under the assumptions of smooth incident fields and p 1 the Neumann traee behaves as 

u^ ~ CjTjfi rj —>• 0, — 1 < Xj < 0, Cj G C where rj denotes the radial distance to the corner Xj 
and Xj is a solution of the transcendental equation 


sin{Xj 7 r — {1 + Xj) 6 j) ^_P + 1 


sin(Aj7r) 


p-1 


(5.4) 


where 9j is the aperture of the interior angle at the corner Xj. In the case of smooth incident 
fields and p = 1, it can be shown that u is aetually more regular. Similar arguments allow 
us to obtain identieal exponents in the corner asymptotie behavior of the solution pN of the SIE 


equations (3.18) 


In what follows we introduce the graded-parameterized version of the four boundary integral op¬ 
erators of the Helmholtz equation. In the light of the discussion above on the regularity of 'Jdu, 
and fi^, we consider the cases when these boundary integral operators act on two types of 2 tt 
periodic densities: (1) we assume that g? G 27r] where a is large enough and in addition (p(t} 

behaves like \t — TjY, r > 0 for all 1 < j < P-|-1; and (2) we assume that G 27r], 0 < /3 < 1 

is a Holder continuous density. We start by dehning the parametrized single layer operator in the 
form 


p 2 n 

{Skg)){t):= Gfc(x(t) - x(r))v9(T)dr. 

Jo 

We define next the parametrized double layer operator in the form 

9Gk{^{t) — x{t)) , 

an(x(r)) 0)\iir)dT. 

and the parametrized weighted adjoint of the double layer operator as 

(K, •m.= f I^WI — ■eO)ir 

Finally, we dehned the parametrized weighted hypersingular operator as 

(' 2 'k 

Jo 

+ PV 




Gfc(x(t) - x(r))|x'(t)| |x'(r)|(n(x(t)) • n(x(T)))V^(T)dr 
[ \x {t)\{dsGk)ix{t) - x(r))V''(r)(iT. 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


Having dehned parametrized weighted versions of the boundary integral operators associated with 
the Helmholtz equation, we present next parametrized weighted versions of the integral equations 


CFIEFK (3.15), CFIESK (3.16), CFIER (3.17), and SCFIE (3.18). 


In the case of the physieal formulations CFIEFK (3.15) and CFIESK (3.16), we introduce the 
additional notation (see B) 

lDu{t) := {^Du){x{t)), := u^^fix{t)), := (77V«*"")(x(t))|x'(t)|, (5.9) 


We multiply both sides of the second equations in formulations CEIEEK ( |3.15[ ) and CFIESK ( 3.16| ) 
by the term |x'(t)| and we obtain 


CFK^ 


7du* 


l.W f 

X := 


-(K 1 +K 2 ) (p-'52 + 5l) 

_-(lVf + pNfi) 







1,W f 

LTaI U1 


^.w rj.inc 

\1n^ \ 


(5.10) 
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and respectively 
CSK“’ 


IDU^ 


( P-' + 1 

I 

+ 

l.W f 

L7a^ ^ j 


1 2 

I 


= 

^w ^.inc 
iN^ 




{K2 - p-^Ki) 


P-HS1-S2) 


l.W f 

L7a1 


(5.11) 


The weighted version of the CFIER formulations (|3.17|) can be written as 

S, 


GFK’' 


IDU^ 

li’i’ t 

L7a1 u 


( P 

-CSK^ + 

2 


\p + i 

l + p 

[-pK 

1 

I 

25J 


P +1 

[-2piV- 

Pl\ 

^,w ^.inc 

UN^ J 


CFK^ 

0 < f < 2 tx. 


IDU- 

L7a1 


(5.12) 


In addition, with efficiency considerations in mind, we also consider an alternative version of the 
regularized formulations in which the operators and are replaced by appropriate Fourier 
multipliers which are the principal symbols of the former operators in the sense of pseudodifferential 
operators [3]. Specifically, we use the following Fonrier multipliers 




nGlj 


and 


{PSs,K4’){xit)) = ^0■5,K(n)^(n)e*”^ crsA^) 






(5.13) 

(5.14) 


acting on 27r-periodic densities (j) and where 4>{n) and '!/>(n) are the Fourier coefficients of the 
functions (j) and "0 respectively. With the aid of these Fourier multipliers we define the weighted 
Principal Symbol CFIER formnlation (CEIERPS) 


PSGFK^ 


7DW 

l.W f 

llN U] 


^ CSK“' + ^ 


p + l 


PSs,. 


l + p [-pPS\ 


N,k 


1 

I 

2 PSs,n] 


p + i 

i-^pPS^,. 

pi \ 



l.W f 

J L 7jv ^ J 


0 < f < 27r. 


(5.15) 


We note that we did not establish the well-posedness of the CFIERPS formulations. Nevertheless, 
we give plenty numerical evidence in Section |5.4| that the CFIERPS formulations are robnst and 
computationally advantageous. 

Einally, in order to derive weighted parametrized versions of the SCFIE equation, we consider the 
same layer representation of the fields and u^, but a weighted Burton-Miller type combination 
of the form for j = I, 2 are 

defined as in (5.2)) to lead to the following weighted boundary integral equation 


where 

and 


^^p^ + - ipS^p^ = - ippDU^' 


K"' = -Rp^ipl - 2K2 ’"’) - pKl ’"'(/ + 2K2 ’"’) + 2 (Arf - N^)S 2 


p£R 7/0, on[0,27r] (5.16) 


= -pSi{I + 2KP'^) -{I- 2Ki)S2. 
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Remark 5.2 It is possible to define more general weighted unknowns if we replaee the weight 
|x'| by Ix'l*^, 5 > 1/2, the new weights leading to new weighted integral equation formulations. 
Furthermore, it is possible to define weighted unknowns that involve the Dirichlet data, that is 
:= ti(x(t))|x'(t)|‘^ leading again to new weighted integral equations. We note that the latter 
type of weighted unknowns (with S = 1/2) and weighted boundary integral equations of CFIESK 
type were used in irnns- 


5.2 Nystrom discretizations based on kernel splitting and trigonometric inter¬ 
polation 


We use a Nystrom discretization of the weighted parametrized equations (5.10), (5.11), (5.12), (5.15), 


and (5.16) that relies on (a) splitting of the kernels of the weighted parametrized operators into 
smooth and singular components, (b) trigonometric interpolation of the unknowns of these in¬ 
tegral equations, and (c) analytical expressions for the integrals of products of periodic singular 
and weakly singular kernels and Fourier harmonics. Several details of this method were originally 
introduced in [28]; for completeness, we give out the full details in what follows. 


We present first a strategy to split the kernels of the weighted parametrized integral operators 


featured in equations (5.10), (5.11), (5.12), and (5.16) into smooth and singular components. The 


latter, in turn, can be expressed themselves as products of known singular kernels and smooth 
kernels. We assume that A: > 0 and we begin by looking at the operator 

r27r /*27r 

{Sk^p){t) = / Mk{t,T)(p{T)dT := / Gfc(x(t) - x(r))v9(T)dr, (5.17) 

Jo Jo 

where (/? it is a sufficiently smooth 27r—periodic function—recall that basically ip{t) := |x'(t)|7)Yu(x(t)). 
The kernel Mk{t,T) can be expressed in the form 


Mk{t,T) = Mk,i{t,T)ln 4 sill 


2 t - T 


-k Mk,2it,T) 


with 


Mk,i{t,T) := - —Jo(fc|r|) 

Mk, 2 {t,T) := Mfc(t,r) - Mfc,i(f,r) In ( 4sin 


2 t - r 


where we have denoted, for lightening this and forthcoming expressions, 

r = r{t, t) = x(t) — x(t). 

Observe that the diagonal terms are given by 


Mk,i{t,t) - 


Mk,2{t,t) = 


_ Jl In 

27r 27r 2 


where C ~ 0.5772 is the Euler-Mascheroni constant. 
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The parametrized double layer operator, see (3.2), is defined as follows 


p27r 

Hk{t,T)'ip{T)dT := 

Jo 


r-27r 




We note that the kernel of the operator behaves as (i) \t — t\~^ when t ^ Tj,t < Tj and r —?■ 
Tj, T > Tj for 2 < j < P and as (ii) (|t —t| mod 27r)“^ when t —>■ Ti = 0 and r —)• Tp+i = 27r (that 
is when x(t) and x(r) approach a corner from different sides). Thus, the integral in the definition 
of the operator Kj. should be understood in the sense of Cauchy Principal Value integral. However, 
for Holder continuous densities ijj, it is possible to recast the operators into an equivalent form 
that features operators that involve integrable expressions only |22) . In order to do so, we express 
Kk in the form 

r2TT 


{Kk'll^){t) = 


+ 


d 


[Gfc(x(t) - x(r)) - Go(x(t) - x(t))] |x'(r)|V^(T)d7 


0 c)n(x(r)) 

clGo(x(t) — x(r)) n \\/ If \ I f w 7 

0 


9n(x(T)) 

Jo 


(9n(x(r)) 


(5.19) 


We note first that 


/o 9n(x(r)) 


1 

2 


if t G [0,27r] \ {Ti,... ,Tp} 


t = Tj,l < j < P, 


where 9j is the inner angle of the jth corner. For the second integral, and since we have assumed V’ 
to be Holder continuous, the integrand is weakly singular. Finally, in the kernel of the first integral 


operator in (5.19) we find the function 

dGk 


Hk{t,T) = 


5n(x(T)) 


ik 




(r) |x'(r)| = jI'(t) -r 


where 


u{t) := n(x(r))|x'(r)| = w' {t){x 2 {w{t)),-x[{w{t))). 
We have in addition the decomposition 

t — t' 


where 


Hk {t, t) = Hk,i [t, r) In 4 sin 


HUt,r) := 

47r r 


+ Hk,2{t,T) 


which have the diagonal terms 

1 v>{t) ■ x"(t) 


Hk,i{t,t) = 0, Hk,2{t,t) = 


4 tt |x'(t)|2 
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It can be easily seen that the second function in the kernel of the first integral operator in (5.19) 
is given by 


Ho{t,T) = 


dGn 


cln(x(r)) 


(r) |x'(t)| = 


1 u{t)■ r 
271 Irp ’ 


Ho{t,t) = 


1 

47 r |x'(t)p ’ 


and thus appearing in the first term in (5.19), can be dehned even at corner 

points where |x'(t)| = 0. 


The graded-parametrized adjoint of the double layer cf. 


is given by 


{Kj, ip){t)= Hk{t,T)(p{T)dT:= |x(t)|-- (p{T)dT. 


dn{x{t)) 


(5.20) 


Here 


ik 


Hk{t,T) = —i'{t)-r 




The kernel [t, r) can be expressed in the form 


{t, t) = Ffc'i(t, r) In ( 4 sin' 


rT / . „■ 2 \ , ttT 




with 




Ji{k\r\) 


and 


■= Hj{t,T)-In (^4 

1 u{t)- x"(t) 




dvr |x'(t)P 


A simple calculation shows that is infinite whenever |x'(t)| = 0, that is w'{t) = 0. 

Remark 5.3 Notice that although H ^2 unbounded in and around corners, the product 
still vanishes at corners. Indeed, in a neighborhood of a corner, say t ~ Tj the function ip behaves 
as ip{t) ~ Cj[w'{t)]^~^^G “1/2 < ^j- ^ careful inspeetion of the singularity o/R/~2(t, t) reveals that 
this expression behaves as w"{t)/w'{t) and thus the product R/'2(t, ~ Cjw"{t)[w'{f)\^G Given 

that w{t) ^ \t — Tj\P, t ^ Tj, we see that the latter product is regular enough for t —)• Tj provided 
that p is sufficiently large. 

Finally, for the graded-parametrized version of the hypersingular operator N^, we add and subtract 
^ ln(4 sin^iit - t)/2) to get 


-j p‘l'K j. _ /*27r /*27r 

cot (t) dr + Qfc(Tr)V’(r)dr-F r)V’'(r) dr 
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with 


Qk{t,r) := fc^Mfc(t,T)(x'(t) • x'(t)) 

Dk{t, t) := ^ In ^sin^ • 


Note we have used 

|x (i)||x'(r)|(n(x(t)) • n(x(r)) = x (t) • x'(r). 

The kernel Qk can be treated similarly to the kernel M^. On the other hand, a simple calculation 
gives that 

Dk{t,T) = Dk,i{t,T)ln ^4sin^ ~T”) 

where 


Dk,2{t^T) 


. r 

Att |r| 

Dkit,T) - Dk^iit,T)ln ^4sin' 


2 t - T 
2 


have diagonal terms 

1 x'(t)-x"(t) 

Dk,i{'t,t) = 0, Dk,2it,t) = — • 

Again, Dk^ 2 {t,t) is infinite at corners, but the trapezoidal rule can still be applied since that term 
is multiplied by ip'it) which vanishes at the corners—this requires the same type of justification 
used in Remark 15.31 


We note that the weighted integral equations CFIESK (5.11) and SCFIE (5.16) feature the differ¬ 
ence operator — NJf ■ While this difference can be performed directly using the methodology 
presented above for the evaluation of the operators Nf and N 2 , a more advantageous approach re¬ 
lies on the methods developed by Kress in [23] for the evaluation of operators — Nq where Nq is 
the weighted hypersingular operator corresponding to wavenumber A; = 0. The latter methodology 
consists of expressing the graded-parametrized operators, constructed from (3.4) instead, as 


r-27r 


iWk - K]'^){y) = - / -Go)(x(t)-x(r))(t,r)i/(T)V’(r)dT (5.21) 

Jo 


We have (see the proof of Lemma 4.2) that 




V^(Gfc-Go)(r) = -L_^W(fc|r|)^rr' 


1 


„T 


r(i)/ 


r - 


27r 


= Lpfc,o(L'r)ln sin 


2 t - r 


+ L 2 ,k,o{'t, T) 


2 T 1 
——r r — T-r7 


where I is the 2x2 identity matrix and 


Ji(A;|r|) 




k 

Li,kAt,r) = ^ 


L2,k,o{^’T) ■= V^(Gfc-Go)(r)-Lpfc,o(L'r)ln sin 


|r| 

t — T 


r r 


23 

















satisfies 


^2 ^ 
Li,k,o{t,t) = —I, L2kfl{t,t) = k 

OTT 

It follows then that 


JL In f 

47r V 2 


i 2C-1 


Svr 

2 t - r 


I + ,Lio x'(t) (x'(t))’^ 


47 r |x'(t)p 


(i/(t))’^(V^(Gfc - Go)(r))iv(r) = Li_fc(t,r)ln ^4sin2 ~T”) 
with diagonal terms 

2C-1 


Li,k{t,t) = —|x'(t)|^ L 2 ,k{t,t) = k'^ 

OTT 


JL In ( 

47r V 2 




Svr 




that are bounded even around corner points where x'(t) = 0. Thus, we can apply the procedure 
above for the graded-parametrized operator = {NJf^ — N^) — — N^) so that we are 

led to integral operators whose kernels are of the form 

Lk^M (t, t) := [Li^ki (t, t) - Lpfcj (t, r)] In ^4 sin^ “T”) ~ ' 

The splitting techniques presented above can be adapted for the evaluation of the operators 5^ 
and N1^ with Q{k) > 0 using additional smooth cutoff function supported in neighborhoods of the 
target points t according to the procedures introduced in jl]. 


5.3 Trigonometric interpolation 


We describe next a Nystrom method based on trigonometric interpolation that follows closely the 
quadrature method introduced by Kress in [23], which in turn relies on the logarithmic quadrature 
methods introduced by Kussmaul [26j and Martensen |29|. The main idea is to use global trigono¬ 
metric interpolation of the quantities that are the solutions of the integral 

equations (5.10), (5.11), and (5.16). Given that the larger the exponent p of the sigmoidal transform 


is, the smoother the quantities and /r"' are, the trigonometric interpolants of these quanti¬ 

ties converge fast with respect to the number of interpolation points. We choose an equi-spaced 
splitting of the interval [0, 27r] into 2n points. We choose Tj such that — Tj are proportional 
(with the same constant of proportionality) to the lengths of the arcs of T from Xj to Xj+i for all 
j. Consequently, the number of discretization points per subinterval [Tj,Tj^i], 1 < j < P may 
differ from each other. We thus consider the equi-spaced collocation points ... ,4nLi} 

and the interpolation problem with respect to these nodal points in the space T„ of trigonometric 
polynomials of the form 


v{t) = 


n 

E' 

m=0 


n—1 


i mt -|- E 


sin mt 


m=l 


which is uniquely solvable [2T|. We denote by Pn '■ C'[0,27r] —)• T„ the corresponding trigonometric 
polynomial interpolation operator . We use the quadrature rules |23| 


c 2 n 


In I 4 sin^ 


t — T 


r-27r 


f{r)dT 


t — T 


/ In ( 4 sin^ 

Jo V 2 

2n-l 

E «!”’(*)/(*!”’) 

i=0 


{Pnf)iT)dT 


(5.22) 
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where the expressions are given by 


0 ?! 1 -j 

R\"'\t) = —- — cos m{t — 1^) —^ cosn(t — 

* n ^ m * n2 * ^ 

m=l 


We also use the trapezoidal rule 

r27r 


r2TT 


>0 


f{r)dT: 


{Pnf){T)dT = fit 


2n-l 




i=0 


Finally, we have the quadrature rule 


1 

— / cot 
dvr In 




/ ir)dT 


1 

— cot 
dvr In 


T — t d 

2 dr 


[iPnf)ir)]dT 


2n-l 






i=0 


where 


n—1 


T^it) = -— E "^cosm(t - 4”^) - -cosn(t - 


m=l 


(5.23) 


(5.2d) 


The derivatives of the densities needed for the evaluation of the operators NJf and are effected 
by differentiation of the global trigonometric interpolant of the densities. This can be pursued 
either by means of Fast Fourier Transforms (FFTs) or using the Fourier differentiation matrix 
whose entries are given by = ^(—1)*+.? cot ^ = 0. 

Finally, given that the values of vf and fjP vanish at corner points, the terms in the boundary 
integral equations that feature these quantities are not collocated at corner points. Alternatively, 
this issue can be bypassed altogether by shifting the mesh by h/2, where h is the meshsize. 
All of the interpolatory quadratures presented above still apply for the shifted meshes. Finally, the 
Fourier multipliers PSg^ and defined in equations (5.13) and (5.Id) can be easily evaluated 

using trigonometric interpolation and FFTs. 


5.4 Numerical results 

We present in this section a variety of numerical results that demonstrate the properties of the var¬ 
ious formulations considered in this text. Solutions of the linear systems arising from the Nystrom 
discretizations of the transmission integral equations described in Section are obtained by means 
of the fully complex, unrestarted version of the iterative solver GMRES [35]. For the case of the 
regularized formulations we present choices of the complex wavenumber k in each of the cases 
considered; our extensive numerical experiments suggest that these values of k lead to nearly op¬ 
timal numbers of GMRES iterations to reach desired (small) GMRES relative residuals. We also 
present in each table the values of the GMRES relative residual tolerances used in the numerical 
experiments. 

We present scattering experiments concerning the following two Lipschitz geometries (see Figure [^: 
(a) a square whose sides equal to d, and (b) a U-shape scatterer of sides equal to d and indentation 
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Figure 1: Domains for our numerical tests: (a) square; (b) U-shape domain 


equal to 2 . For every scattering experiment we consider plane-wave incidence of direction 
(0,-1) and we present maximum far-field errors, that is we choose sufficiently many directions 
(1024 directions have been used in our computations) and for each direction we compute the 
far-field amplitude defined as 


ni(x) 



«oo(x) +0 



X —oo. 


(5.25) 


The maximum far-held errors were evaluated through comparisons of the numerical solutions 
corresponding to either formulation with reference solutions by means of the relation 


Soo = max |u^“^^(i) - {k)\ (5.26) 

The latter solutions were produced using solutions corresponding with rehned discretizations 

based on the formulation CFIESK with GMRES residuals of 10“^^ for all geometries. Besides far 
held errors, we display the numbers of iterations required by the GMRES solver to reach relative 
residuals that are specihed in each case. We used in the numerical experiments discretizations 
ranging from 6 to 12 discretization points per wavelength, for frequencies ki and k 2 in the medium 
to the high-frequency range corresponding to scattering problems of sizes ranging from 2.5 to 40.8 
wavelengths. The columns “Unknowns” in all Tables display the numbers of unknowns used in each 
case, which equal to the value 4n dehned in Section for the weighted GEIEEK, GEIESK, CEIER, 
and GEIERPS formulations, and 2n for the weighted SGEIE formulation. In order to remind the 
reader of the fact that the SGEIE formulations require half the number of unknowns required by 
each of the other formulations, we denote the former by SGEIE* in the tables. Eollowing common 
practice [3|, we used the GEIEEK operators as their own preconditioners, and we denote these 
by GEIEEK^. We note that in this case the computational time required to perform a matrix- 
vector product corresponding to the GEIEEK^ formulation is double that related to the GEIEEK 
formulation. 


We start by presenting the high-order convergence of our Nystrom solvers in Tables [T]^ We have 
used sigmoid transforms with a value p = 3 for all formulations but SGEIE, in which case we used 
p = 4. The need for a different value of p for the latter formulations can be attributed to the more 
singular nature of the solutions of these equations. As it can be seen, solvers based on the GEIESK 
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ki 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 


Iter. 

^OO 

Iter. 


1 

4 

256 

32 

6.6 X 10^^ 

34 

2.4 X 10^“ 

43 

3.8 X 10^“ 

43 

3.2 X 10^^ 

32 

3.1 X 10^® 

1 

4 

512 

31 

8.0 X 10^^ 

34 

1.8 X 10^'^ 

46 

5.5 X lO^'* 

43 

3.9 X lO-'* 

33 

3.7 X lO-'* 

1 

4 

1024 

31 

1.0 X 10^^ 

34 

1.1 X 10^® 

49 

8.0 X 10^® 

47 

4.8 X 10^® 

34 

4.6 X 10^® 

1 

4 

2048 

31 

1.2 X lO^^' 

34 

4.1 X 10^1° 

54 

1.2 X 10^® 

47 

6.0 X 10^® 

34 

5.8 X 10^® 


Table 1: Scattering experiments for the square geometry with p = 1, and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected rj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k2)(2 + i ki. The GMRES 
residual was set to equal 10“^^. 


fei 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^CO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

1 

4 

352 

84 

4.3 X 10-^ 

75 

3.7 X 10-^ 

64 

1.4 X 10^^ 

73 

1.7 X 10^^ 

62 

1.9 X 10^^ 

1 

4 

704 

82 

5.2 X 10^® 

75 

2.2 X 10-® 

66 

1.7 X 10^® 

74 

2.1 X 10^® 

63 

2.3 X 10-® 

1 

4 

1408 

81 

6.4 X lO-'* 

75 

1.5 X 10^® 

67 

2.4 X lO-'* 

75 

2.6 X lO-'* 

63 

2.9 X lO-'* 

1 

4 

2816 

80 

7.9 X 10^® 

75 

9.9 X 10^® 

68 

3.7 X 10-® 

77 

3.1 X 10^® 

63 

3.7 X 10-® 


Table 2: Scattering experiments for the U-shape geometry with p = 1, and for the CFIEFK^, 
CPIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected rj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k^)!"! + i ki. The GMRES 
residual was set to equal 10“^^. 


formulations are the most accurate on account of the facts that they do not require numerical 
differentiation. 

In Table we present computational times required by a matrix-vector product for each of the for¬ 
mulations CEIEEK, CFIESK, SCFIE, CFIER, and CFIERPS. The computational times presented 
were delivered by a MATLAB implementation of the Nystrom discretization on a MacBookPro 
machine with 2 x 2.3 GHz Quad-core Intel i7 with 16 GB of memory. We present computational 
times for the square geometry, as the computational times required by the U-shaped geometry 
considered in this text are extremely similar to those for the square geometry at the same levels 
of discretization. As it can be seen from the results in Table the computational times required 
by a matrix-vector product for the CFIEFK, CFIESK, SCFIE, and CEIERPS formulations are 
quite similar, while the computational times required by a matrix-vector product related to the 


fei 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^CO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

1 

4 

256 

58 

9.9 X 10-'‘ 

39 

1.5 X 10^® 

48 

2.1 X 10-® 

60 

2.0 X lO-'^ 

76 

4.1 X 10-'‘ 

1 

4 

512 

56 

1.2 X lO-"* 

39 

9.0 X 10-'® 

49 

3.3 X 10--* 

52 

4.5 X 10-® 

80 

5.2 X 10-® 

1 

4 

1024 

54 

1.5 X 10^® 

37 

6.0 X 10-® 

51 

5.0 X 10-® 

57 

6.0 X 10-® 

84 

6.5 X 10-® 

1 

4 

2048 

53 

1.9 X 10^® 

37 

4.1 X 10-** 

52 

7.6 X 10-® 

57 

7.0 X 10-'^ 

87 

8.2 X 10-'^ 


Table 3: Scattering experiments for the square geometry with p = and for the CFIEFK^, 

CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected p = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k^^jl -|- i ki. The GMRES 
residual was set to equal 10“^^. 


27 















































h 

k2 

Unkno-wns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

1 

4 

352 

no 

6.5 X 10^'“ 

67 

4.8 X 10^^ 

71 

5.4 X 10-5 

93 

3.5 X lO-'* 

115 

2.5 X lO-'* 

1 

4 

704 

107 

1.0 X lO-"* 

64 

1.1 X 10^3 

71 

8.0 X 10-"* 

86 

7.2 X 10-5 

119 

3.4 X 10-5 

1 

4 

1408 

107 

2.0 X lO-'^ 

64 

2.5 X lO-'^ 

72 

1.2 X 10-* 

88 

1.3 X 10-5 

123 

8.1 X 10-5 

1 

4 

2816 

105 

3.9 X 10-® 

63 

5.7 X 10-5 

72 

1.7 X 10-5 

91 

4.0 X 10-5 

126 

3.4 X 10-5 


Table 4: Scattering experiments for the U-shape geometry with p = and for the CFIEFK^, 

CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected p = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k2)(2 + i ki. The GMRES 
residual was set to equal 10“^^. 


Unknowns 

CFIEFK 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

256 

4.5 

4.8 

4.2 

5.6 

5.0 

512 

15.9 

16.4 

15.4 

19.1 

17.0 

1024 

59.4 

63.6 

64.2 

73.0 

63.8 


Table 5: Computational times (in seconds) for the matrix-vector products (seconds) needed by 
the formulations CEIEEK, CFIESK, SCEIE, CFIER, and CFIERPS for the square geometry. 


CFIER formulation are on average 1.16 times more expensive than those required by the other 
three formulations. 

We present next in Tables |6[|^ the performance of our solvers based on the five formulations con¬ 
sidered in this text in the case of high-contrast, high-frequency configurations. We conclude, 
in conjunction with the results presented in Table that solvers based on the SCFIE, CEIER, 
and CFIERPS formulations consistently outperform solvers based on the classical formulations 
CFIEFK^ and CEIESK in the regime under consideration. Furthermore, the solvers based on the 
SCFIE and CFIER formulations compare favorably to the solvers based on the CFIESK formula¬ 
tion. 


We conclude with an illustration in Tables 10 -11 of high-contrast high-frequency scenarios whereby 
the computational gains associated with solvers based on the SCFIE, CFIER, and CFIERPS are 
the most significant. As it can be seen in Table 11, solvers based on the SCFIE formulations and 
CFIER formulations can result in computational times that are at least 3 times faster than those 
based on the classical CFIESK and CFIEFK formulations. 


ki 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

1 

4 

256 

18 

6.8 X 10-5 

24 

3.7 X 10-* 

25 

3.8 X 10-5 

26 

3.1 X 10-5 

21 

3.1 X 10-5 

2 

8 

512 

24 

2.9 X 10-5 

39 

5.0 X 10-* 

37 

2.1 X 10-* 

33 

3.6 X 10-5 

32 

1.4 X 10-5 

4 

16 

1024 

62 

4.7 X 10-5 

94 

4.0 X 10-5 

63 

1.0 X 10-* 

58 

7.5 X 10-5 

62 

1.8 X 10-5 

8 

32 

2048 

119 

8.1 X 10-5 

162 

8.2 X 10-5 

112 

3.8 X 10-* 

102 

6.6 X 10-5 

115 

6.7 X 10-5 


Table 6: Scattering experiments for the square geometry with p = 1, and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected p = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k2)/2 + i ki. The GMRES 
residual was set to equal 10“^. 
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ki 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCEIE* 

CEIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^OO 

1 

4 

352 

52 

4.3 X 10-2 

64 

7.0 X 10-3 

45 

1.4 X 10-2 

49 

2.0 X 10-2 

44 

2.2 X 10-2 

2 

8 

704 

76 

2.5 X 10-2 

107 

3.2 X 10-3 

78 

9.4 X 10-3 

79 

1.5 X 10-2 

75 

2.0 X 10-2 

4 

16 

1408 

117 

8.7 X 10-3 

149 

7.7 X 10-3 

136 

3.3 X 10-3 

124 

3.9 X 10-3 

113 

4.4 X 10-3 

8 

32 

2816 

281 

3.3 X 10-2 

351 

3.2 X 10-2 

257 

2.8 X 10-3 

257 

1.4 X 10-2 

244 

1.2 X 10-2 


Table 7: Scattering experiments for the U-shape geometry with p = 1, and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected rj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k 2 )/2 -\- i ki. The GMRES 
residual was set to equal 10“^. 


ki 

k2 

Unknowns 

CFIEFK2 

CFIESK 

SCEIE* 

CEIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^oo 

Iter. 

Coo 

Iter. 

Coo 

1 

4 

256 

30 

9.8 X 10-3 

23 

6.7 X 10-3 

28 

2.1 X 10-3 

33 

2.1 X 10-3 

47 

4.3 X 10-3 

2 

8 

512 

53 

1.2 X 10-3 

34 

1.3 X 10-3 

59 

3.6 X 10-3 

39 

6.5 X 10-3 

72 

8.5 X 10-3 

4 

16 

1024 

82 

6.7 X 10-3 

53 

7.0 X 10-3 

88 

2.5 X 10-3 

51 

6.7 X 10-3 

99 

1.8 X 10-3 

8 

32 

2048 

236 

1.6 X 10-3 

112 

2.1 X 10-3 

205 

2.3 X 10-3 

111 

1.8 X 10-3 

197 

4.3 X 10-3 


Table 8; Scattering experiments for the square geometry with p = A:f//c 2 , and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected tj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki-\- k 2 )/2 -\- i ki. The GMRES 
residual was set to equal 10“^. 


ki 

k2 

Unknowns 

CFIEFK2 

CFIESK 

SCEIE* 

CEIER 

CFIERPS 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

1 

4 

352 

76 

5.4 X 10-3 

40 

4.8 X 10-3 

57 

5.5 X 10-3 

58 

5.1 X 10-3 

66 

3.5 X 10-3 

2 

8 

704 

98 

7.7 X 10-3 

64 

2.0 X 10-3 

87 

9.6 X 10-3 

66 

1.9 X 10-3 

107 

8.1 X 10-3 

4 

16 

1408 

236 

1.7 X 10-3 

126 

2.2 X 10-3 

168 

4.3 X 10-3 

128 

9.2 X 10-3 

181 

2.5 X 10-3 

8 

32 

2816 

424 

3.2 X 10-3 

252 

2.8 X 10-3 

286 

7.0 X 10-3 

216 

1.1 X 10-3 

305 

3.3 X 10-3 


Table 9: Scattering experiments for the U-shape geometry with p = kf/k'^, and for the CFIEFK^, 
GFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected tj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki-\- k 2 )/2 -\- i ki. The GMRES 
residual was set to equal 10“^. 


ki 

k2 

Unknowns 

CPIEFK3 

CFIESK 

SCEIE* 

CFIER 

CEIERPS 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

Iter. 

Coo 

3.5 

1 

256 

16 

2.1 X 10-3 

23 

8.1 X 10-3 

21 

1.1 X 10-3 

20 

1.7 X 10-3 

21 

3.3 X 10-3 

7 

2 

512 

29 

1.7 X 10-3 

41 

2.0 X 10-3 

24 

4.4 X 10-3 

20 

2.1 X 10-3 

30 

2.0 X 10-3 

14 

4 

1024 

59 

8.9 X 10-3 

56 

1.7 X 10-3 

35 

3.2 X 10-3 

22 

1.1 X 10-3 

57 

4.2 X 10-3 

28 

8 

2048 

85 

4.1 X 10-3 

94 

9.4 X 10-3 

39 

5.5 X 10-3 

25 

1.1 X 10-3 

87 

1.2 X 10-3 


Table 10: Scattering experiments for the square geometry with p = 1, and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected p = ki. In 
the regularized formulations CFIER and CFIERPS we used k = {ki + k2)/2 + i 4. The GMRES 
residual was set to equal 10“^. 
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fel 

k2 

Unknowns 

CFIEFK^ 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^00 

Iter. 


Iter. 

^00 

Iter. 

^00 

Iter. 

^00 

3.5 

1 

352 

42 

3.4 X 10-^ 

51 

2.4 X 10-^ 

30 

3.6 X 10-^ 

26 

1.5 X 10-^ 

37 

2.4 X 10-3 

7 

2 

704 

54 

3.2 X 10^® 

70 

2.7 X 10-2 

35 

3.0 X 10-1 

33 

1.3 X 10-3 

47 

2.7 X 10-3 

14 

4 

1408 

123 

2.7 X 10-2 

148 

9.5 X 10-2 

47 

3.2 X 10-1 

46 

1.2 X 10-3 

77 

2.6 X 10-3 

28 

8 

2816 

238 

9.8 X 10-2 

240 

1.3 X 10-1 

84 

3.6 X 10-1 

67 

1.8 X 10-3 

169 

4.2 X 10-3 


Table 11: Scattering experiments for the U-shape geometry with p = 1, and for the CFIEFK^, 
CFIESK, SCFIE, and CFIERPS formulations In the SCFIE formulation we selected rj = ki. In 
the regularized formulations CFIER and CFIERPS we used k = (fci + k2)/2 + z 4. The GMRES 
residual was set to equal 10“^. 


ki 

k2 

Unknowns 

CFIEFK2 

CFIESK 

SCFIE* 

CFIER 

CFIERPS 

Iter. 

^OO 

Iter. 

^OO 

Iter. 

^00 

Iter. 

^00 

Iter. 

^00 

1 

4 

256 

20 

5.0 X 10-3 

25 

2.7 X 10-1 

19 

8.7 X 10-3 

30 

5.8 X 10-3 

21 

5.8 X 10-3 

2 

8 

512 

24 

2.8 X 10-3 

41 

3.7 X 10-1 

30 

1.4 X 10-3 

32 

5.6 X 10-3 

32 

1.7 X 10-3 

4 

16 

1024 

62 

4.0 X 10-3 

97 

1.8 X 10-3 

55 

1.4 X 10-3 

59 

5.6 X 10-3 

62 

1.8 X 10-3 

8 

32 

2048 

122 

7.9 X 10-3 

173 

5.6 X 10-3 

96 

3.3 X 10-3 

103 

8.7 X 10-3 

117 

4.6 X 10-3 


Table 12: Scattering experiments for the Bq, q = 512 sphere of radius 2 with p = 1, and for the 
CPIEFK^, CFIESK, SCFIE, and CEIERPS formulations In the SCFIE formulation we selected 
r] = ki- In the regularized formulations CFIER and CFIERPS we used n = {ki + k2)/2 + i ki. The 
GMRES residual was set to equal 10“^. 


Finally, we present in Table [T^ a comparison between scattering solutions corresponding to Lipschitz 
scatterers and solutions corresponding to nearby smooth scatterers that are obtained from rounding 
the corners. More specifically, we considered the sphere of radius 2 in using the norm for 
q = 512, that is 


Bn 


Jq := {(xi, X 2 ) G M : xf + x| = 2'^}, q = 512 

which is a close and smooth (rounded) approximation of the square geometry in Figure [I] —indeed, 
the distance between the scatterers is about 10“^. We compare the far-field signature of Bq,q = 
512 to that of the square for various wavenumbers using the various boundary integral equation 
formulations considered in this text. We note that for a given frequency the numbers of iterations 
required by boundary integral formulations to reach the same tolerance are roughly the same for 
the square and its very close smooth (rounded) approximation—see Tables]^ and 12 


6 Conclusions 


In this work we have presented high-order Nystrom discretizations based on polynomially graded 
meshes for several boundary integral formulations including certain regularized formulations for 
Helmholtz transmission problems in domains with corners. We have rigorously proven the well- 
posedness of some of these formulations and have shown that solvers based on the regularized and 
on the single integral equations outperform solvers based on commonly used boundary integral 
equation formulations in the case of high-contrast, high-frequency applications. The numerical 
analysis of these schemes will be subject of future investigation. Extensions of the regularization 
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scheme used in this paper to the case of multiple dielectric scatterers will also be subject of future 
investigation. 
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